First-Order Superconducting Transition of Sr 2 Ru0 4 
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By means of the magnetocaloric effect, we examine the nature of the superconducting-normal (S-N) transition 
of Sr 2 Ru04, a most promising candidate for a spin-triplet superconductor. We provide thermodynamic evidence 
that the S-N transition of this oxide is of first order below approximately 0.8 K and only for magnetic field 
directions very close to the conducting plane, in clear contrast to the ordinary type-II superconductors exhibiting 
second-order S-N transitions. The entropy release across the transition at 0.2 K is 10% of the normal-state 
entropy. Our result urges an introduction of a new mechanism to break superconductivity by magnetic field. 



The order of a phase transition provides one of the most 
fundamental pieces of information of the long-range ordered 
state accompanied by the phase transition. In case of super- 
conductivity, the order of the superconducting-normal (S-N) 
transition in magnetic fields reflects how the superconductiv- 
ity interacts with the magnetic field and how it is destabilised. 
For example, for a type-I superconductor, the in-field S-N 
transition is a first-order transition (FOT) [1], because of an 
abrupt disappearance of the superconducting (SC) order pa- 
rameter caused by the excess energy for magnetic-flux exclu- 
sion. For a type-II superconductor, in contrast, the in-field S-N 
transition is ordinarily a second-order transition (SOT) [1], In 
this case, penetration of quantized vortices with accompany- 
ing kinetic energy due to orbital currents leads to a continuous 
suppression of the SC order parameter up to the upper critical 
field H C 2- This type of pair-breaking is called the orbital effect. 

A well-known exception for type-II superconductivity is the 
case where the superconductivity is destroyed by the Zeeman 
spin splitting [2]. When the spin susceptibility in the SC state, 
Xsc, is lower than that in the normal state, Xn, the SC state 
acquires higher energy AEz ~ (l/2)(^ n - Xsc)noH 2 with re- 
spect to the normal state, due to the difference of polarizabil- 
ity of the electron spin. This destroys superconductivity at the 
Pauli limiting field /d Hp ~ [2/ioEcond/Ctn -;tsc)] 1/2 , where 
AEz reaches the SC condensation energy E con d. Such a pair- 
breaking effect is called the Pauli effect. It is theoretically 
predicted that a strong Pauli effect leads to a first-order S-N 
transition at temperatures sufficiently lower than the critical 
temperature T c [3]. This prediction has been confirmed in a 
few spin-singlet superconductors [4—6]. 

The type-II superconductor S^RuC^ (T c = 1.5 K) is one 
of the most promising candidates for spin-triplet supercon- 
ductors [7-9]. Due to its unconventional superconducting 
phenomena originating from the orbital and spin degrees of 
freedom as well as from non-trivial topological aspect of the 
SC wave function, this oxide continues to attract substantial 
attention [10-14]. The spin-triplet state has been directly 
confirmed by extensive spin susceptibility measurements by 
means of the nuclear magnetic resonance (NMR) using sev- 
eral atomic sites [15-17] and the polarized neutron scatter- 
ing [18]: Both experiments have revealed ^- sc = Xn in the 
entire temperature-field region investigated. This means that 
Hp oc (Xn ~Xac) ^ 2 is infinite and the Pauli effect is irrelevant 
in this material. 



Interestingly, several properties of the S-N transition of 
Sr2RuC<4 have not been understood for more than 10 years 
within the existing scenarios for the spin-triplet pairing. For 
example, H C 2(T) is more suppressed than the expected behav- 
ior for the orbital effect, when the field is parallel to the con- 
ducting ab plane [19-21]. In addition, several quantities such 
as the specific heat C [20], thermal conductivity k [20], mag- 
netization M [22], exhibit sudden recovery to the normal-state 
values near H C 2 for H \\ ab and at low temperatures. 

To resolve the origin of such unusual behavior, we per- 
formed measurements of the magnetocaloric effect (MCE) of 
Sr2RuC<4. The MCE is a change of the sample temperature T 
in response to a variation of the external magnetic field H; we 
measure T while sweeping H at a constant rate. The thermal 
equation of the MCE is written as [23] 
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where S is the entropy, C is the heat capacity of the sample, 
k is the thermal conductance between the sample and thermal 
bath, H is the sweep rate of the magnetic field, T\,- d ± is the 
temperature of the thermal bath, and d'gioss is the dissipative 
loss of the sample. When k is small so that the second term 
is negligible, the equation reduces to the relation for the con- 
ventional adiabatic MCE. In the other limit where the thermal 
coupling between the sample and bath is strong, the first term 
in turn becomes negligible, leading to the "strong-coupling 
limit" relation [23] (dS/dH) T = -(kAt/H) - T-\d'Q ioss /dH) 
with At = (T - 7b a th)/r [24]. In this limit, the measured 
At is linearly dependent on (dS/dH)r- Thus, it is expected 
that T and At exhibit peak-like anomalies at a FOT and step- 
like anomalies at a SOT. Because of this qualitative difference, 
the strong-coupling MCE is suitable to distinguish a FOT and 
a SOT. We found that our calorimeter indeed works nearly 
in this strong-coupling limit, with the first term in Eq. (1) 
amounting to at most 10% of the second term. We however 
didn't neglect the first term in the evaluation of the entropy 
discussed below. 

For the present study, we used single crystals of Sr 2 Ru04 
grown by the floating-zone method [25]: Sample #1 weighing 
0.684 mg with T c = 1 .45 K and Sample #2 weighing 0. 1 84 mg 
with T c = 1 .50 K. The value of T c of Sample #2 is equal to the 
ideal T c of Sr2Ru04 in the clean limit [26], indicating its ex- 
treme cleanness. The MCE was measured using a hand-made 
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FIG. 1. (color online) (a,b) Representative raw data of the mag- 
netocaloric effect (MCE) of Sr 2 Ru0 4 for T ~ 0.2 K at H || ab 
(H ~ [100]). The dotted curves indicate t^To for the up sweep 
(pink) and the down sweep (cyan), which corresponds to the back- 
ground contribution [24]. (c,d) Relative temperature change —Atf 
(red) and At± (blue) due to the MCE for the same field condition at 
different temperatures. For clarity, each pair of curves is shifted ver- 
tically by 0.25% and 0.4% for panels (c) and (d), respectively. The 
high-temperature (T > 0.6 K) data in (d) are multiplied by 5. 



sensitive calorimeter. Magnetic field was applied using a vec- 
tor magnet system [27]. Details of the experimental method is 
described in the Supplemental Material [24]. 

We first present the MCE for H || ab (H ~ [100]) and 
T ~ 0.2 K measured at fioH = +1.02 mT/sec in Figs. 1(a) 
and (b). Obviously, T(H) exhibits peak-like behavior near 
H C 2, rather than a single step-like behavior. This feature 
becomes clearer in the background-subtracted At^{H) (up- 
sweep) and Ati(H) (down-sweep) curves shown in Figs. 1(c) 
and (d) [24]. The observed peak provides indication of a FOT 
in S^RuCv Note that a slight asymmetry in the MCE sig- 
nal (i.e. \At-^(H)\ < \Ati(H)\) is attributed to the energy dissi- 
pation mainly due to vortex motion causing a heating in both 
the field up-sweep and down-sweep measurements [28]. More 
importantly, H C 2 is clearly different between the up-sweep and 
down-sweep curves. The difference between the up-sweep 
onset // C 2f and the down-sweep onset H C 2[ is approximately 
H{)AH C 2 = /uo(H C 2i - H C 2[) = 20 mT for Sample #1 and 15 mT 
for Sample #2. This difference corresponds to 15-20 sec for 
fioH = 1.02 mT/sec. The difference cannot be attributed to an 
extrinsic delay of the temperature measurement, since the de- 
lay time of our apparatus is much shorter than 15-20 sec [29]. 
We have also confirmed that a finite AH C 2 is observed for 
lower sweep rates such as fioH = +0.2 mT/sec. Therefore, 
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FIG. 2. (color online) Field angle 6 variation of the magnetocaloric 
effect of Sample #2 at T ~ 0.2 K. The red and blue curves indicate 
—Atf(H) andAti(H), respectively. Each curve is shifted vertically by 
0.5% and the data for 6 > 2.0° are multiplied by 3 for clarity. 



this difference of H C 2 is indeed intrinsic, and provides defini- 
tive evidence that the S-N transition is a FOT accompanied by 
supercooling (or possibly superheating). Note that the very 
sharp peak in At{H) at H C 2 for Sample #2 demonstrates the 
cleanness and homogeneity of this sample. 

Next, we focus on the variation of the MCE with temper- 
ature and field angle. As represented in Figs. 1(c) and (d), 
both the peak in At(H) and the supercooling becomes less pro- 
nounced as temperature increases. Around 0.8 K, these fea- 
tures totally disappear and the S-N transition becomes a SOT 
as expected for ordinary type-II superconductors. In Fig. 2, 
we present several MCE curves for fields tilted away from the 
ab plane toward the c axis by the amount which we define as 9. 
When the field is tilted only by ~ 2 degrees, the FOT features 
disappear. 

From the MCE data for H || ab, we deduce the entropy us- 
ing Eq. (1) [24]. Figure 3(a) again characterizes the FOT with 
a huge peak in (dS ldH)j and supercooling/superheating. In 
Fig. 3(b), we present AS = S - S a = $"(dS /8H) T dH di- 
vided by temperature. Here, S n is the entropy in the normal 
state. The total entropy S can be calculated with the assump- 
tion SJT = y e , where y e = 37.5 mJ/K 2 mol is the electronic 
specific heat coefficient [30]. The jump in S/T across the FOT 
is approximately 6S /T = -3.5 ± 1 mJ/K 2 mol at the low- 
est measured temperatures. This value of 5S /T amounts to 
approximately 10% of S n /T, and the latent heat L = TSS 
at 0.2 K is 0.14 + 0.04 mJ/mol. We can check the con- 
sistency of this value using the Clausius-Clapeyron equation 
yUod// C 2/dr = -5S/5M, where 5M is the jump in M across 
the FOT. Using the values fi dH c2 /dT ~ -0.20 ± 0.05 T/K 
estimated from our H c %\ data for Sample #2 and 5M ~ 
-0.014 emu/g from the magnetization study [31], we obtain 
6S/T - -5.2 + 1.2mJ/K 2 molfor0.2K. This value reasonably 
agrees with the value from our MCE experiment. In addition, 
S at lower fields also exhibits agreement with other thermo- 
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FIG. 4. (color online) Superconducting phase diagram of Sr2Ru04 
for H || ab (H ~ [100]) deduced from the magnetocaloric effect for 
Sample #2. The red squares and the blue crosses indicate the onset 
H C 2 for the up- and down-sweeps, respectively. The inset presents 
H c2 and AH cl (triangles) in the low-temperature region. 



FIG. 3. (color online) (a) Field dependence of (jxqT) ~'(dS /dH) T of 
Sr 2 Ru04 deduced from the magnetocaloric effect, (b) Field depen- 
dence of AS IT . In both (a) and (b), the main panels present data for 
Sample #2 and the insets for Sample #1; the solid and broken curves 
present up- and down-sweep data, respectively. The right axis of (b) 
indicates S IT obtained by assuming S n /T = 37.5 mJ/K 2 mol [30]. 
The double-headed arrow in (b) illustrates the jump 6S/T = -3.5 + 
1 mJ/K 2 mol at the transition. 



dynamic studies [20, 22], as explained in the Supplemental 
Material [24]. 

We summarize the present observations in the phase dia- 
grams presented in Figs. 4 and 5. The region for which the 
FOT emerges is limited to temperatures below Tfot ~ 0.8 K 
for 6 = 0° and field angles within \6\ < 2° for T ~ 0.2 K. 
Interestingly, the FOT region is included in a wider region in 
which the behavior of H Q i cannot be described solely by the 
conventional orbital effect [21]: H C 2(T) substantially deviates 
from the linear behavior and H C 2(&) cannot be fitted with the 
effective mass model (Fig. 5). These facts indicate that the 
ordinary orbital effect cannot be a origin of the FOT. 

Let us compare the present results with previous observa- 
tions. The rapid recoveries of k/T [20], and M [22] near H C 2 
for H || ab have been observed in the region where the S- 
N transition is revealed to be of first order. Thus, it now turns 
out that these recoveries are actually consequences of the FOT. 
However, supercooling (or superheating) at the S-N transition 
in Sr2Ru04 has never been reported in previous studies. This 
is probably because the supercooled metastable normal state 
easily nucleates into the SC state. Thus, a fast and continu- 
ous sweep is helpful to observe the supercooling, rather than 
point-by-point measurements. The smallness and cleanness 
of the present samples have also assisted the observation, be- 
cause the number of nucleation centers (e.g. surface defects, 
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FIG. 5. (color online) Field-angle dependence of (squares) and 
H c2 i (circles) at T ~ 0.2 K. The green curve indicates the fitting to 
the data of Sample #1 in the range \8\ > 2.0° [21] with the effective 
mass model H c2 (6) = // c2 (90°)/(sin 2 6 + cos 2 8/T 2 ) 1 ' 2 , where T is the 
anisotropy parameter and is obtained to be 25 from the fitting. 



lattice imperfections) is reduced for small and clean samples. 
In contrast to the previous studies on the bulk SC phase, a hys- 
teresis in the in-field S-N transition was observed for the in- 
terfacial 3-K phase superconductivity in the Sr2RuO/t-Ru eu- 
tectic [32, 33]. Possible relation between this hysteresis and 
the present observation is worth examining. We note that the 
second transition revealed by the C/T measurement [20] at 
H2, which is 20-30 mT below H C 2, cannot be attributed to the 
onset of the FOT. Although we have not so far obtained con- 
vincing MCE data supporting the H2 anomaly, we need higher 
experimental resolution to clarify this issue. 

In the rest of this Letter, we discuss the origin of the 
FOT. As we have explained, the FOT should originate from 
a pair-breaking effect beyond the conventional orbital effect. 
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Naively, a possible candidate of such a pair-breaking effect 
is the Pauli effect [34]. However, in the case of Sr2Ru04, 
NMR and neutron studies have revealed Xsc — Xn [15-18]. In 
particular, the observed negative hyperfine coupling provides 
strong evidence that NMR correctly detects the spin suscepti- 
bility [9]. In addition, x%c = Xn nas been confirmed for differ- 
ent nuclei at several atomic sites. Other experiments also indi- 
rectly support the spin-triplet scenario [10-12, 14]. Therefore, 
the Pauli effect should be absent in Sr2Ru04 and the FOT can- 
not be attributed to the Pauli effect either. Thus, an unknown 
pair-breaking effect, or in other words, a non-trivial interac- 
tion between superconductivity and magnetic field, must be 
taken into account. 

For a spin-triplet superconductor with Xsc = Xn, it is naively 
expected that the \H) condensate can be mutually converted to 
the Iff) condensate by magnetic field without destroying the 
SC state. However, in contrast to the expectation, the present 
experiment shows that the triplet SC state in S^RuC^ is no 
more stable at high fields where the Zeeman splitting is no 
longer a perturbation. Indeed, H C 2 for T — » nearly matches 
the field Hp = {2poE COD d/xsc) 1 ^ 2 ~ 1.4 T, where the Zeeman 
spin energy in the SC state (1 /2)x S cP()H 2 is equal to E con & [35]. 
This fact suggests that the Zeeman splitting between |f T) an d 
111) condensates, which has not been considered in the exist- 
ing theories on the SC phase diagram of Sr2Ru04 [34, 36-39], 
is a source of the non-trivial coupling between magnetic field 
and triplet superconductivity. 

Let us propose possible mechanisms of the non-trivial in- 
teraction. We can categorise them into microscopic and 
macro/mesoscopic mechanisms. The microscopic mecha- 
nisms include the pinning of the electron spin direction 
at certain Appoints predicted by the band calculation [40] 
and confirmed by the angle-resolved photoemission spec- 
troscopy [41]. Such a pinning of the spin direction may lead to 
a constraint in the spin-polarization due to the Zeeman effect. 
The closeness of the Fermi energy to the van Hove singular- 
ity [42] is also worth considering, because a slight modifica- 
tion of the chemical potential due to the Zeeman effect might 
disturb the pairing glue. 

The macro/mesoscopic mechanisms include possible inter- 
actions among the Cooper-pair orbital angular momentum L, 
the Cooper-pair spin S, the vortex vorticity, and the magnetic 
field. Indeed, a pair-breaking effect due to L was proposed 
in Ref. [43], although this theory cannot be directly applied 
as long as the orbital motion is assumed to be purely two 
dimensional. As another macro/mesoscopic mechanism, the 
kinematic polarization discussed in the context of the stability 
of the half-quantum vortex (HQV) is instructive [14, 44]. It 
was proposed that a velocity mismatch between |fT) and |XX) 
condensates around a HQV results in a shift of the chemical 
potential of these two condensates due to difference in their 
kinetic energies and leads to an additional spin polarization 
coupling to the magnetic field. By an analogy to this theory, 
we expect that consideration of kinematics of the condensates 
in high fields may provide a route to unveil the non-trivial 
coupling between Cooper-pair and magnetic field. 



In summary, our MCE study of S^RuC^ revealed defini- 
tive evidence for a first-order S-N transition in the low- 
temperature region for fields nearly parallel to the ab plane. 
The FOT, not attributable to conventional mechanisms, indi- 
cates a non-trivial interaction between spin-triplet supercon- 
ductivity and magnetic field. This new information on the 
bulk superconductivity serves as a basis for investigations of 
the non-trivial topological nature of the SC wavefunction as- 
sociated with the "Majorana-like" edge modes. We also antic- 
ipate that the abrupt growth of the order parameter across the 
FOT, accompanied by vortex formation and non-trivial sym- 
metry breaking, should provide a new playground for inves- 
tigation of novel vortex dynamics, which might be related to 
quantum turbulence and/or to the Kibble-Zurek mechanism. 
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Shingo Yonezawa, Tomohiro Kajikawa, Yoshiteru Maeno 
Department of Physics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan 

I. DIFFERENCE BETWEEN THE ADIABATIC AND STRONG-COUPLING MAGNETOCALORIC EFFECTS 



As already mentioned in the main text, the basic formula for the 
magnetocaloric effect (MCE) is expressed as: 







\dHj T 





AT 
AH 



k(T - r bath ) 

TH 



1 d'Qioss 
T AH ' 



(SI) 



External field H 
(Sweep rate H) 

thermometer 1 



sample 



Entropy S 
Heat capacity C 
Temperature T 



electrical wire 



Thermal conductance k 



bath 



Temperature T 



where S is the entropy, C is the heat capacity of the sample, k is 
the thermal conductance between the sample and thermal bath, H 
is the sweep rate of the magnetic field, T is the sample tempera- 
ture, Tbath is the temperature of the thermal bath, and d'gioss is the 
dissipative loss of the sample (see Fig. SI). We call the case where 
the second term is negligibly small as the adiabatic MCE (AMCE) 
and where the first term is negligibly small as the strong-coupling 
MCE (SMCE). 

To demonstrate the difference in the behavior of T(H) between 
the AMCE and SMCE, we solve Eq. (SI) numerically by assuming 

functional forms of S (H), C(H), k(H) and A'Q\ otifi (H). For simplicity, we here only treat the case k(H) = const, and A'Q\ osti (H) = 0. 
We further assume C(H) = const, in order to extract the essential difference between the AMCE and SMCE, although C in reality 
should depends on H as well as on T, especially near a phase transition. For S (H), we adopt simple linear field dependences to 
simulate a first-order transition (FOT) and a second-order transition (SOT) as shown in the top four panels of Fig. S2: 



FIG. SI. Schematic drawing of a typical setting for magne- 
tocaloric effect measurements. 



S fot(H) - 



A(H - H c ,) + B(H C * 
B(H-H C ) + S {) , 
So, 



H c ) + S, 



(H < H„) 
(H ct <H<H C ) 
(H c < H); 



(S2) 



and 



S sm(H) 



A(H-H C ) + S , 
So, 



(H < H c ) 
(H > H c ). 



(S3) 



In these equations, H c denotes the critical field of the phase transition and H c * is the onset field of the FOT introduced to simulate 
a realistic imperfection with a finite broadening of the FOT. Here, we used the FOT width H c - H ct = 0.01H C . The coefficients 
A and B are the slopes of S(H) in the low-field phase and within the FOT region, respectively. We here assume B/A = 10, i.e. a 
10-times steeper slope within the FOT region than in the low-field phase. The entropy in the high-field phase is assumed to be a 
constant So- 

The results of the calculation for different values of k are presented in Fig. S2. As the sample-bath relaxation time C/k 
becomes smaller, the system evolves from the AMCE to SMCE. As one can see, T(H) is almost proportional to (dS ldH)j for 
the SMCE as expected from Eq. (SI). Thus, the qualitative difference between a FOT (peak in T(H)) and a SOT (step in T(H)) 
is quite clear, although the MCE signal T(H)/ Tbath is relatively small. In contrast, for the AMCE, qualitative difference between 
the shapes of a FOT curve and of a SOT curve becomes less pronounced, whereas the available MCE signal can be very large 
compared to that of the SMCE. 

To conclude this section, we demonstrate that the SMCE is a reliable method to distinguish a FOT from a SOT in spite of 
relatively small changes in the sample temperature. For the SMCE, such a determination can be done just from the qualitative 
shape of the raw T(H) curve even without detailed analyses of the curve. 
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FIG. S2. Demonstration of differences between the adiabatic magnetocaloric effect (AMCE) and the strong-coupling magnetocaloric effect 
(SMCE) for a first-order transition (left column) and a second-order transition (right column). The assumed functional form of 5 (H), as well as 
its field derivative (dS /dH) T , is presented in the top panels. The calculated MCE responses with different values of the sample-bath relaxation 
time C/k are presented in the rest of the panels. The red and blue curves indicate field-up sweeps and field-down sweeps, respectively. The 
system approaches the AMCE limit for larger C/k, whereas it approaches the SMCE limit for smaller C/k. The parameters H and C are chosen 
to be H/H c = 0.001 sec"' and CT bMh /(AH c ) = 1.0. Note that, for C/k = 10 3 sec, a deviation from the adiabatic limit shows up, for example, 
in the field dependence of T for H > H c . 



II. DETAILS OF THE EXPERIMENTAL METHOD 



For the present study, we used single crystals of Sr2Ru04 grown by the floating-zone method [SI]: Sample #1 weighing 
0.684 mg with T c = 1.45 K and Sample #2 weighing 0.184 mg with T c = 1.50 K. The value of T c of Sample #2 is equal to the 
ideal T c of Sr2Ru04 in the clean limit [S2], indicating its extreme cleanness. The samples were cut, cleaved, and polished; their 
T c was checked by heat capacity and AC susceptibility measurements. We found that Sample #1 have mosaic structure with a 
c-axis tilting of ~ 0.2°, while Sample #2 is free from such mosaicity. The mosaicity in Sample #1 does not affect our conclusion. 

We developed a sensitive calorimeter for MCE measurements consisting of a small thermometer and a heater made of com- 
mercial thick-film RuC<2 resistors fixed with low thermal conductance Pt-W wires. To avoid complication due to the SC transition 
of solder, we removed the solder coating of the resistors and used silver paste to attach the wires to the resistors. The calorimeter 
was cooled with a dilution refrigerator. We measured k and C separately using the relaxation and ac methods. Magnetic field was 
applied using a vector magnet system consisting of two orthogonal magnets and a horizontal rotating stage [S3]. For the MCE 
measurement, we mounted the crystal with its c-axis nearly horizontal. We used the horizontal magnet only because the system 
does not allow simultaneous field change of the two magnets. Thus, the accuracy of the field alignment is better than 0.1° with 
respect to the ab plane but is ~ 5-10° for the azimuthal direction within the ab plane. In order to improve the signal-to-noise ratio 
of the MCE curves, we repeated field sweeps for 20-30 times and take averages among up and down sweep curves separately. 
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III. PROCESS OF THE EVALUATION OF THE ENTROPY 

In this section, we describe details of the process to evaluate the entropy S(H) from the MCE data shown in Fig. 1 of the main 
paper. 



Equation (SI) can be rewritten as 



dS_ 

dH 



C (AAt\ kAt 



(S4) 



where At = (T - T^^/T and Di oss (T, H) — (l/T)(d'Qi oss /dH). We need to evaluate At(H) and Di oss from the raw data in order 
to obtain the entropy S (H). 
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FIG. S3. Description of the process of the entropy evaluation at T ~ 0.21 K and fioli = 1.02 mT/sec for Sample #1. For all data shown here, 
the red curves present up-sweep data and blue curves present down-sweep data, (a) Raw data of T and T bi (up sweep: magenta, down sweep: 
cyan). Tq = T(1.45 T) and 7b g o = 7b g (1.45 T) are indicated with the crosses, (b) Normalized temperatures T/Tq and 7bg/7bgo. (c) Reduced 
temperature At = [(T/T ) - (7 , bs /r bs0 )](r /T ). (d) Heat capacity of the sample plus addenda measured by the relaxation method. The broken 
curve presents the result of the fitting. This fitted values were used for the entropy-evaluation process, (e) Nominal entropy derivatives f[ and 
£l containing the contribution from the loss term D\ osi . The broken curves indicate the contribution of the first term (i.e. the AMCE term) of 
the right-hand side of Eq. (SI). The fact that this contribution is at most 10% of and indicates the present calorimeter nearly works in the 
strong-coupling limit, (f) Loss term D loss = (£j - ff)/2. The broken curve is the interpolation into the FOT region, (g) Field derivative of the 
entropy (dS/dH) T . (h) Entropy evaluated by integrating (dS/dH) T . 



In the evaluation process of At(H), we managed to solve an issue on reading errors of thermometers. A temperature reading 
r read , represented in Fig. S3(a), inevitably contains a certain reading error 5T with respect to the true temperature r tme : i.e. 
r read = r true + 6T. The error 6T may originate from errors in the calibration function ((?r cal ) and from long-term (~ days) 
drifts in the reading in the electronic devices (6T dntt ). In the present case, since 5T can be on the same order as the intrinsic 
temperature variation due to the MCE, we need to eliminate errors originating from 6T for accurate evaluation of S (H). 
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The only way to cancel out unknown calibration errors dT cal in T and 7b a th is to use a temperature measured with the same 
thermometer as T, instead of 7b at h- Thus, we adopt 7b g , which is the sample temperature for fields away from the ab plane. In 
the present study, 7b g was measured in fields 20° away from the ab plane; such a tilt of the magnetic field reduces H C 2 down to 
0.2 T at 0.2 K. Using 7b g , the extrinsic (i.e. background and normal-state) contribution of the MCE signal can be subtracted and 
we can evaluate the entropy change only due to the superconductivity. 

To minimize contributions from 6T dnft , we re-define At using the normalized temperatures T/To and 7b g /7bg0 (Fig. S3(b)) as 

T(H)/T -T b JH)/T b „ 

At(H) = -J-JL^ bgV bgU , (S5) 

T(H) T 



with To = T(H = Hq) and 7b g o = T bs (H = Ho), where Ho is a certain field H = Ho above H C 2 for H || ab. In the case of Fig. S3, 

^true 



we chose fXoHo = 1.45 T. We have indeed confirmed that the reading errors are canceled out and At ^ 1 - 77 me /7 ,true U p to the 



first order in small values, as we explain in detail in Appendix. The obtained At is presented in Fig. S3(c). 

Other two important quantities, C and k, are determined from separate measurements using the relaxation-time method and 
the ac method. The obtained C data (Fig. S3(d)) is consistent with the previous works [S4, S5, S6]. We found that k was field- 
independent within our experimental resolution in the present field range and we adopt k — 1 1.4 nW/K for the analysis of the 
data in Fig. S3. 

Next, we need to evaluate the loss term Di 0S s> which leads to an asymmetric (\Atf\ + \At§) MCE signal. For the present case, 
dealing with a type-II superconductor, Di oss originates from heating due to incoming and escaping motion of vortices associated 
with both increasing and decreasing fields. For the present case, the MCE signal is indeed asymmetric (i.e. \At^\ < |z)fj) in the 
superconducting (SC) state, even down to fields much smaller than H C 2 (see Fig. 1, and Fig. S3). Note that, in general, a FOT 
is often accompanied with an energy dissipation [S7]. However, for the present case, it is difficult to conclude that the energy 
dissipation characteristic of the FOT has any substantial contribution to the MCE signal in this case, because we did not detect 
additional asymmetry near the FOT as shown in Fig. S3(f) 

The loss term Di oss is estimated by taking an average of the up-sweep and down-sweep results: We first evaluate the "nominal" 
entropy derivatives 

?TB -_£_(^)-*£ (S6) 
tT \-Aty\dHj H 



for up-sweep data (H > 0) and 



C (dAt,\ kAh 

^-i^Mrir (S7) 

for down-sweep data (H < 0), as shown in Fig. S3(e). Then Di oss (Fig. S3(f)) is obtained by 

Ao SS (7 , , J ff) = i(^-f T ). (S8) 

This process is quite similar to that used in Ref. [S8]. This is based on the assumption that Di oss is independent of the sweep 
direction. We should be careful, here, that the right-hand side of Eq. (S8) can be finite even if Di oss = near a FOT because of 
supercooling/superheating. Thus, we need to estimate Di oss within the FOT regions by an interpolation using a simple polynomial 
function as indicated by the broken curve Fig. S3(f). This interpolation is performed so that the entropy-conservation law is 
satisfied between above and far below H C 2- 

With the estimated A oss , (OS /8H)t can be obtained using Eq. (S4) as shown in Fig. S3(g). Then, as presented in Fig. S3(h), 
the entropy S (H) can be calculated by the integration 



s{H) - s ^Cl^l m (s9) 



where S n is the normal-state entropy. Note that the evaluations of (dS /dH)r and S become less accurate when Di oss dominates 
the MCE. 



IV. COMPARISON OF SAMPLE #1 AND SAMPLE #2 



In order to demonstrate the validity of our entropy evaluation, we here compare results for two Sr2Ru04 crystals described in 
the main text: Sample #1 (the SC transition temperature T c = 1.45 K with a broader transition in magnetic fields) and Sample #2 
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(T c = 1.50 K with an extremely sharp transition even in magnetic fields). When the loss term Di oss is smaller than the intrinsic 
contributions, we can accurately evaluate (dS ldH)j as well as S both for Sample #1 and Sample #2. As presented in Figs. S4(a) 
and (c), (dS ldH)j of both samples agrees with each other below 1 .25 T at 0.4 K and below 1 .05 T at 0.8 K, although deviations 
originating from the difference in H C 2 and the sharpness of the SC transition have been observed near H C 2- Accordingly, the S (H) 
curves, drawn with an assumption that the normal-state entropy S n J = 37.5 mJ/K 2 [S9] is common between Sample #1 and #2, 
reflect the similarity and difference in (dS ldH)j- The data at 0.3 K shown in Fig. S5 also exhibits good agreement considering 
the characteristics of each sample. This agreement between different samples indicates the validity of our analyses. 
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FIG. S4. Comparison of the entropy S(H) and its field derivative (dS/dH) T for Sample #1 and #2, at (a,b) 0.4 K and (c,d) 0.8 K. In all panels, 
the red and blue curves indicate the values for Sample #1, and the magenta and cyan curves indicate the values for Sample #2. The S(H) 
curves in (b) and (d) are obtained with an assumption that the normal state entropy 5„ is 37.5 mJ/K 2 mol [S9] for both samples, as indicated 
with the dotted lines. 

At lower temperatures, the contribution of £V SS becomes more significant. For Sample #2, since Di oss dominates the MCE 
signals as shown in Fig. 1 of the Main Text, it is quite difficult to accurately evaluate the intrinsic (dS /dH)r except near H C 2, 
where the intrinsic contribution is still larger than A oss . Indeed, the jump in the entropy across H C 2, 6S /T ~ 3.5 mJ/K 2 mol, is 
consistent between the two samples, as already mentioned in the Main Text. 



V. COMPARISON WITH OTHER THERMODYNAMIC STUDIES 



In this section, we compare the entropy obtained from our MCE results, Smce, with those from other thermodynamic studies. 
The specific heat C has a relation to the associated entropy Sc as 




(S10) 



To evaluate Sc, which are plotted in Fig. S5(a), we first extrapolated the C(T)/T data in Ref. [S4] to K using a polynomial 
function. We then integrated C/T up to 0.32 K. We also plot (dSc/dH)r in Fig. S5(b) obtained by simply taking two-point slopes 
of the Sc(H) data, together with {dSy[CEldH)j- As demonstrated in Figs. S5(a) and (b), Sc and Smce, as we H as (dSc/dH)r and 
(dSMCE/dH)r, reasonably agrees with each other in the present field range. 
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FIG. S5. Comparison of the MCE results with previous thermodynamic studies, (a) Comparison of S(H) obtained from the MCE and specific 
heat [S4]. (b) Comparison of (dS/8H) T obtained from the MCE, magnetization [S10], and specific heat [S4]. 



The Maxwell's relation yields the relation between the magnetization M and the associated entropy Sm as 



dS M \ (dM\ 

mj T = ^[dT) H 



(SH) 



Therefore, we can estimate {dSMldH)j from M reported by Tenya et al. [S10]. Because they reported M(H) curves only at a 
few fixed temperatures, we estimated (dSM/dH)r as follows: We first extracted M(T) data from Fig. 2 of Ref. [S10], and then 
fitted a quadratic function M(T) = Mo +AT 2 to the data. This fitting was successful only below LIT; and note that, even below 
1.1 T, only 3-5 data points are available for the fitting. Because {dMjdT)u ~ 2AT, (dSM/dH)r divided by ^lqT is given by 2A, 
which is plotted in Fig. S5(b) after an appropriate conversion of the unit. In spite of the uncertainty due to the small number of 
data for the fitting, (8SM/dH) T and (_dSucEldH) T reasonably agree with each other. 



VI. SUMMARY 

To summarize this Supplementary Material, we explain our careful entropy evaluation process in detail and demonstrate that 
the process is valid although the T variation due to the MCE is less than 1%. We emphasize that a large MCE signal (AT(H)) 
may not necessarily mean that the evaluated entropy is accurate. As temperature variation due to the MCE becomes larger, a 
violation of the constant-temperature condition becomes much serious. Then, for a substantially large MCE signal, evaluation of 
S(H) would become complicated because it leads to a large error if one evaluates S(H) by integrating (dS /8H)t from a single 
measurement curve. In contrast, in our approach, we can safely integrate the experimentally-obtained (dS /8H)t to obtain S (H) 
because the constant-temperature conditions is almost satisfied. Therefore, our approach provides an alternative route to evaluate 
the entropy from MCE results. 



Appendix: Details of the estimation of reading errors in thermometry 

Readings of the measured sample temperature r lead inevitably contains reading errors ST with respect to the true temperature 

y^true. 



(SA1) 
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This is also the case for Tq, T bg and 7bgo- We here evaluate influence of these errors to the modified definition At — [(T/Tq) - 
(T bg /T bg o)]/(T /To). Putting Eq. (SA1) to the definition of At (Eq. (S5)), the experimental value of At can be written as 

At=l-^-% (SA2) 

7 bgO 1 

T^ + 6T hgT «™ + 6T 



T hoa + ST b°o T true + 6T 



bgO 
fttae Ttme 

I - bg 
Ttrue Ttrue 
bgO 



ST bg ST bg0 ST ST 
1 + — rr^ : — + 



ytrue Ttine Ttrue Ttrue 
bg bgO J 



(SA4) 



Here, we neglect second order terms of the small values. Because of the approximation 

ST 5T ST L T tme -T Q me \ ST 



1 - 



Time Ttrue , /Ttrue _ Ttrue\ T-true I Ttrue I Ttrue 

1 '0 > \ I 1 Q 

within the first order to the small values (such as ST/T me and (r true - T Q me )/T me ), and because of similar approximations for 
7b a and 7b g o, we can substitute the denominators in the parentheses in Eq. (SA4) by T^ as 

r bT 77™ / 6T bg - rJTbgo -ST + ST \ 
At^l--£ — l + — ^ -. (SA5) 

Ttrue Ttrue 1 Ttrue I v ' 

bgO * ' 

The reading errors ST can be separated into a term originating from calibration errors and a term originating from electronics 
drift: 

ST = ST C *\T, H) + ST drift . (SA6) 

Note that we can use the same function ST cd[ (T, H) for all four temperatures discussed here because they are measured with the 
same thermometer. In addition, the drift term should satisfy ST dntt = ST^ ft and 8T^ n = ST^. Recalling the fact that Tq and 
Tbgo are the values at a certain field H - Hq, we can write the reading errors as 

ST = ST c * l (T tme , H) + ST dm , (SA7) 
ST = ST c *\T me , H ) + (5r drift , (SA8) 
<Sr bg = ST cal (T^ e , H) + ST d f, (SA9) 

r-drif 
bg 



ST bgf) = ST™\T l ™, Ho) + ST d f. (SA10) 



The temperature variation in the calibration error <5r caI should be small in the present small temperature range, i.e. 
ST CR \T tme , H) = ST CR \T l ™ e ,H) and ST cal (T me , H ) =* ST cd \T l ™, H ). Corrections to this simplification only results in higher 
order errors. Thus, we can conclude 

ST bg - ST bg0 - ST + ST (S A 1 1 ) 

and 

Time j, trae Ttrue 

At « 1 - — 2_ = 1 - J2L (v r?™ = TSf) (SA12) 

Time Ttme Ttrue v U bgu 7 v 7 

bgO J J 

up to the first order in small values. 

To summarize, use of the modified definition/If = [(T/To) - {T bg /T bg o)] / (T /To) minimizes influences of the reading errors to 
the entropy evaluation. 
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